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of the hehcal spin crystal: a candidate for the partially ordered state of MnSi 
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MnSi is an itinerant magnet which at low temperatures develops a helical spin density wave. 
Under pressure it undergoes a transition into an unusual partially ordered state whose nature is 
debated. Here we propose that the helical spin crystal (the magnetic analog of a solid) is a useful 
starting point to understand partial order in MnSi. We consider different helical spin crystals and 
determine conditions under which they may be energetically favored. The most promising candidate 
has bcc structure and is reminiscent of the blue phase of liquid-crystals in that it has line-nodes of 
magnetization protected by symmetry. We introduce a Landau theory to study the properties of 
these states, in particular the effect of crystal anisotropy, magnetic field and disorder. These results 
compare favorably with existing data on MnSi from neutron scattering and magnetic field studies. 
Future experiments to test this scenario are also proposed. 
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MnSi is a well studied itinerant ferromagnet T] but 
has recently been at the focus of renewed attention since 
the discovery of unusual properties in high pressure stud- 
ies 1^ H, Q, 1^. The ambient-pressure magnetic phase 
of MnSi is characterized by a hierarchy of three major 
energy scales 6:|. First are interactions that favor itin- 
erant ferromagnetism. At a much lower energy scale, 
Dzyaloshinskii-Moriya (D-M) spin-orbit coupling (which 
is allowed in the non-centrosymmetric B20 crystal struc- 
ture of MnSi) leads to a spiraling of the magnetic mo- 
ment with a single helicity (helical spin-density wave). 
The small ratio between these scales is responsible for 
the long spiral pitch (A = 180A vs. the lattice constant 
a = 4.6^). The D-M interaction defines chirality and 
the length of the ordering wave- vector {Q — 2tt/\) but 
not its direction. The latter is pinned to (111) by even 
weaker crystal-anisotropy terms [(111): class of crystal 
directions related by cubic symmetry]. Elastic neutron 
scattering at atmospheric pressure reveals Bragg spots 
at these eight points in the Brillouin zone which are at- 
tributed to domains each containing a single spiral state. 
On application of pressure however, a first order phase 
transition occurs at pc = 14.6kbar. At higher pressures, 
unusual properties are observed. First, the neutron scat- 
tering signal is completely changed from the low pressure 
phase. Enhanced scattering is seen at wavevectors with 
a length similar to the low pressure phase, but with in- 
tensities no longer sharply peaked along (111). Rather, 
the intensity is more diffuse over the wavevector sphere, 
hence the name "partial order", but clearly peaking at 
(110) 121 . Secondly, the resistivity displays a non fermi 
liquid (NFL) temperature dependence 4] . Although they 
both onset at Pc, the NFL persists to pressures far be- 
yond where "partial order" is seen. The relation between 
these two puzzles is therefore unclear - here we focus only 
on offering a theory for the "partially ordered" state. Re- 
cent theoretical work on collective modes and electronic 
properties of helimagnets are in Refs. [tI. Other the- 
oretical proposals for the high pressure state of MnSi 



have invoked proximity to a quantum multi-critical point 
8], magnetic liquid-gas transitions 9] and skyrmion-like 
structures [l(lllll|. 

None of these theories has attempted to explain the 
new peak positions. As a first attempt, one might spec- 
ulate that the crystal anisotropy which pins the spiral 
along (111) at low pressure is modified at high pressure 
and pins it along a different set of directions, (110). How- 
ever, this view is untenable on many counts, e.g. it is at 
odds with the magnetic field studies 0- Moreover, the 
usual crystal anisotropy terms consistent with the lattice 
symmetry give energy minima either at (111) or (100); 
but (110) are always saddle points (which can be proved 
more generally |l2l|). Similarly, while it is tempting to 
interpret the partially ordered phase of MnSi entirely in 
terms of a dircctionally fluctuating spin spirals, this re- 
quires both a novel theoretical framework for its descrip- 
tion, as well as an explanation for the unusual anisotropy. 

In this work, we study ordered magnetic states which 
are linear superpositions of spin spirals with different 
wave-vectors. We show that in addition to the three en- 
ergy scales (ferromagnetic, D-M and crystal anisotropy) 
that are usually considered in MnSi, there is an addi- 
tional scale that is important, the interaction between 
different modes (expected to be intermediate between the 
D-M and crystal anisotropy), which determines the spin 
crystal state that is selected. We propose that the tran- 
sition under pressure between the single-spiral and par- 
tially ordered states is driven by a change in the inter- 
mode interactions that goes from preferring a single- 
spiral ground state to a multi-spiral ground state. 

Magnetic weak crystallization theory: The Landau free 
energy to second order in the local spin density M(r) of 
a system with full rotation symmetry (transforming both 
space and spin together) - but no inversion symmetry is 



F2 = (roM^ -f J(VM)2 -f 2D M • (V X M)) 



(1) 



where (...) means averaging over the sample and ro, J, D 
are parameters (J > 0). Coupling to fermions which 
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leads to damped dynamics is ignored in the following 
since we will focus on ground states with long range or- 
der. The last term of Eq. is the D-M interaction. 
Clearly, the energy is minimal for circularly polarized 
waves of fixed helicity which satisfy D m* • (iq x niq) = 
— jfj g |mqp, where mq are the Fourier modes of M(r). 
Equation then becomes F2 = J2q''^{'l) l™qP) where 
''(9) = '"0 — D'^/J + J{q — |_D|/J)^. Here, in analogy to 
standard liquids (and in contrast to the uniform ferro- 
magnetic case), r(g) is minimized at a finite wavevector 
Q — \D\/J and for r{Q) 0, all modes on the surface 
|q| = Q in reciprocal space become soft (as opposed to a 
single point q = 0). To study the implications of such a 
singular surface on the phase transition itself is beyond 
the scope of this letter. But we recognize that when 
r{Q) 0, the Gaussian theory not only becomes unsta- 
ble towards formation of a helical spin-density wave along 
any direction, but that linear combinations of such states 
(helical spin crystals) are equally natural candidates. 

In the spirit of weak crystallization theory |l3l| , we now 
study minima of the free energy in the ordered phase, 
where r(Q) < 0. The degeneracy between a simple spin 
spiral and linear combinations of several spin spirals (spin 
crystals) is lifted by interactions i.e. by the fourth-order 
term (F4) in M (odd terms in M are forbidden by time 
reversal symmetry). We assume that F4, as F2, has full 
rotation symmetry and we will include the weak crystal 
anisotropy last. Hence, with q4 — — (qi + q2 + qs): 

F4 = ^ [/(qi, q2, qa) (niq^ • niqj (niqj • niqj . (2) 

qi.q2,q3 

If the interaction is strictly local in real space [F4 cx 
(M^), which is equivalent to a constant coupling 
t^(qii 12, qa) — Uo], it is easy to show that the single- 
mode spin spiral is the absolute minimum of F2 + F4. In 
the spin spiral state, is constant and it is therefore 
possible to minimize both F2 and F4 independently. 

However, f/(qi,q2,q3) is generally not constant. As- 
suming F4 is small, in that its main effect is to provide an 
interaction between the modes which are degenerate un- 
der F2, the important terms of F4 are those with |qi| = 
1^2! = Iqal = |q4| = Q- Then the couphng function U de- 
pends only on two independent angles 29 = arcc os( qi -q^) 
and 1^/2 = arccos[(q2 - qi) • qa /(I - qi • q2)] This 
mapping allows and to be interpreted as the polar and 
azimuthal angles of a sphere and the coupling U{d, cj)) is 
a smooth function on that sphere satisfying: U{d, cj)) = 
U{tt — 9,(f)) = U{9,2tt — (j)). It is natural to expand U 
in spherical harmonics (Yim) which satisfy this relation 
and wc just retain terms with I < 2. Thus U{9,(f>) = 
Uo + Uii sin 9cos(t) + U20 (3 cos^ 6* - 1) f/22 sin^ 9 cos 2(f>. 

To find the absolute minimum of F2 + F4 for a gen- 
eral angle-dependent interaction would be a formidable 
task. Instead, we first limit ourselves to linear combina- 
tions of six spin spirals, corresponding to wavevectors 
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FIG. 1: Mean-field phase diagram of the six-mode model for 
r{Q) < as a function of the interaction parameters for U20 = 
and Uo > 0. In the gray region, F4 < and higher-order 
terms should be added for stability. 



±kj (for j = 1,...,6), where ki = Q/V2 (1, -1, 0), 
k2 = Q/V^ (-1,-1,0), etc. (see inset of Fig. [TJ. We 
may write the six modes as rxikj = Tna*_^ = ^pj ^ where 

nq = [qx (zxq)+izxq]/[l — (z-q)^]^/^. The choice of the 
unit vector z is arbitrary. A different choice leads to a set 
of phase changes in the six complex variables ipj . The in- 
teraction Eq.[21 written in terms of the i/'- variables, con- 
tains terms of the form Vjj' IV'iPlV'i' P (three parameters) 
and one term A 3? {T^ +Ty + Tz), where = '4>i4'2'4'b'4'6, 
Ty = V'i'02V'3V'4 ^iid Tz = —'03 ■04 "05 "06 which is sensitive 
to relative phases. These are the only quartic terms in- 
variant under the microscopic symmetries of the prob- 
lem: translations, time reversal and tetrahedral point 
group operations. The parameters Vjj' and A can be eas- 
ily expressed as a linear combination of J7o, U20 and 
U22- The resulting mean-field phase diagram is shown in 
Fig. n Different ordered phases are separated by first- 
order phase boundaries. In the "spiral" state, only one 
out of six amplitude is non-zero. The square lattice state 
"□" is an superposition of two orthogonal spin spirals 
with equal amplitude. In the region "A" , there are actu- 
ally two degenerate states, one with 1-0^1 > (i = 1, 3, 5) 
and one with (i = 2, 4, 6) (remaining amplitudes be- 
ing zero in each case). Finally in BCCl (A > 0) and 
BCC2 (A < 0), all six modes contribute with equal am- 
plitudes. These states have the periodicity of a BCC 
crystal and will be the focus of this paper. Note that a 
small correction of the purely local interaction C/q, given 
by Uii « U22 ~ 0.2J7o is sufficient to induce a transition 
from the spiral state to the BCCl spin crystal, [lij 

Properties of the BCC spin crystal: From the sole as- 
sumption (suggested by experiments) that a spin texture 
has equal- weighted Bragg peaks at the set of wave- vectors 
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±kj, weak crystallization theory restricts an infinite va- 
riety of physically very different magnetization patterns 
M(r) down to only two possibilities: BCCl and BCC2. 
In BCCl (BCC2), the relative phases of tpi, ■ ■ ■ ,tpe are 
locked in such a way that Tx,Ty and are all negative 
(positive), i.e. for BCCl: ^pi — ±iip4ip5, ip2 — T^4V'6' 
"03 = V'4'/'5V'6 where ipj = ipj/lipjl and for BCC2: ipi = 
±i'04-05, -02 = zLiipi^l^Q, ^1^3 = —'ipi'ipB'ipQ. Three phases 
(e.g. "04, ■05, -06) are arbitrary due to global translation 
symmetry. In addition to the translational degeneracy, 
both BCC states are two-fold degenerate due to time- 
reversal symmetry breaking. In contrast to the single 
spiral, time-reversal [M(r) — M(r)] is not equivalent 
to a translation in either of the BCC states. In fact, both 
BCC states feature a macroscopic time-reversal symme- 
try breaking order parameter S = {M^MyMz) ^ 0. 

BCCl can be defined by its symmetry properties as 
being the unique structure that is invariant under time 
reversal (sign change of M) followed by a 7r/2 rotation 
about the x, y or z axis. For BCC2, the same operations 
result in translations. The real space BCCl structure is 
illustrated in Fig. [3 There are two sets of straight lines, 
along which M is constrained by the symmetry proper- 
ties of BCCl. First, the magnetization vanishes along the 
X, y and z axes [and their translations by (^, ^, ^)], due 
to the invariance mentioned above. These are remines- 
cent of the blue phases of chiral nematic liquid crystals 
[l^ . where line defects also exist but are protected by 
topology, rather than by symmetry. Second, in the cen- 
ter between four parallel node lines, the magnetization 
direction is constrained as shown in Fig. 13 (a). Since 
these properties are determined by symmetry, they are 
stable, even if higher Fourier modes are included. In 
contrast, BCC2 has no node lines, but six point nodes in 
each primitive unit cell. Fig. |2 (c) shows the distribution 
of the local magnetization |M| over the sample. 

Anisotropy and locking of crystal directions: So far, 
our model free energy [Eqs. and jSJ] has been com- 
pletely rotation invariant. By choosing those six modes 
which lead to a BCC crystal, we assume that full rota- 
tion symmetry is spontaneously broken, but a global ro- 
tation still leaves the energy invariant. This degeneracy 
is broken by small anisotropic terms. For the MnSi B20 
crystal structure, the leading anisotropic term in powers 
of M and q is of the form: Fa — a^- g{'kj)\'tpj\'^ where 

f;(k) = + ky + k*. The function g has its maxima and 
minima along (100) and (111), respectively, while (110) 
is a saddle point. It is therefore impossible that this 
term locks a single-spiral to (110). From the observed 
spiral orientation along (111) at low pressures, we obtain 
a > 0. Assuming there is no sign change with pressure, 
we find that this locks the magnetic BCC crystal (BCCl 
or BCC2) to the atomic B20 crystal in such a way that 
all kj point along (110). Hence it naturally explains the 
neutron scattering peaks at the locations seen in Ref. ^. 
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FIG. 2: (Color online) (a) Magnetization pattern M(r) of 
BCCl. The black lines are nodes, where M = 0. Along red 
lines, the magnetization direction is as indicated. The blue 
dashed line indicates the location of the cut shown in (b) 
where vectors denote the in-plane magnetization. The nodes 
lines are the centers of anti-vortices and the directed red lines 
are the centers of meron configurations, (c) Probability dis- 
tributions ((5(Af — |M(r)|)) of three different magnetic states: 
a single-spiral state of amplitude Mo and both BCC states 
with IV'il^ = Mq/G. Landau theory predicts this 1/6 ratio of 
amplitudes at the phase boundary between these states. A 
finite-width (5-function was used for the numerical evaluation. 



Magnetic field: A uniform magnetic field couples lin- 
early to the q = mode of the magnetization, m = (M) . 
Eq. ||2J) couples this additional mode to the spin crys- 
tal modes and produces terms of the form m^|0jp, 
(m • kj)^|0jp and one term /^h^ • m, where is cu- 
bic in the variables ■ The behavior of a single spin 
spiral in a magnetic field is characterized by a strongly 
anisotropic susceptibility, which effectively orients the 
spiral axis along the field Once the spiral is ori- 

ented, the susceptibility is comparatively large. In the 
BCC spin crystal states, there is no such anisotropy in 
the linear response, because symmetry does not allow for 
an anisotropic susceptibility tensor. Since the BCC spi- 
rals, unlike the single spiral, cannot be optimally oriented 
along the field, one expects the following sequence of 
phases as a function of increasing magnetic fields: "BCC" 
"spiral" "spin polarized" , in good comparison with 
experiment This expectation is true within our the- 
ory if ^ is not too large. In fact, the term /i -m induces 
a change of the relative amplitudes and phases of the six 
interfering spirals as a function of the magnetic field and 
therefore adds to the ability of the BCC state to adjust 
to an external field. For example, a magnetic field in z- 
dircction suppresses |0i | and enhances |02 1 relative to the 
other four amplitudes [or the opposite, depending on the 
signs of /i and the time-reversal index 5' = {M^MyMz)]- 
This effect should be observable by neutron scattering 
in a single-domain sample, which may be obtained via 
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proper field-cooling. In practice this needs some care 
since the energy splitting between the two time reversed 
states is only cubic in the field: oc H^HyH^. 

Magneto-transport: The broken time-reversal symme- 
try in the BCC spin crystals, but with the absence of a 
uniform magnetization, can give rise to unusual magneto- 
transport in single-domain samples. The symmetric part 
of the conductivity tensor is allowed a linear field depen- 
dence due to time reversal symmetry breaking [S ^ 0): 



Cab = C^O'^ab + aS \tabc\Hc + O (i?^) 



(3) 



E.g. for a field applied along z, the conductivity along 
(x±y) /v^ would display Zmear magneto-resistance. Sim- 
ilarly, the Hall conductivity (the antisymmetric part of 
cr) is allowed an unusual quadratic contribution: 

^ab - <y"eabcH, + a'Seabc\€cde\HdH, + O {H^) . (4) 

Thus, the Hall effect will in general not simply switch 
sign if the direction of the magnetic field is inverted. 

Effect of Disorder: Although the available MnSi crys- 
tals are very clean from the electrical resistivity point of 
view, helical magnetic structures are sensitive to disorder 
at a much longer length scale. Hence disorder effects need 
to be studied. Non-magnetic disorder which couples to 
the magnetization squared, Fdis — J dr Vrfii (r)|M(r)p, 
will not affect single-spiral states which have a constant 
magnitude of magnetization, but they do affect spin crys- 
tal states which have a modulated magnitude [see Fig|3 
(c)]. Therefore the neutron scattering signal of the spin- 
crystal state is expected to have more diffuse scattering 
than the single spiral, consistent with the experimental 
observation that the high pressure phase has diffuse scat- 
tering peaked about (110) while the low pressure phase 
has sharper spots. Anisotropic spreading of the Bragg 
spots parallel and perpendicular to the wavevector sphere 
is also anticipated (13] . Disorder in D=3 is expected to 
destroy true long range order of the spin crystal, leading 
to either a Bragg glass or a disordered state. In either 
case, time reversal symmetry breaking of BCC 1 remains, 
implying that there must be a finite temperature phase 
transition on cooling into this phase. 

NMR and fiSR: The spatially modulated magnetiza- 
tion magnitude of spin-crystal states [Fig. [21 (c)] should 
be visible in muon spin rotation (/xSR) and zero field 
NMR experiments. While /iSR on MnSi has been per- 
formed only at low pressures zero field NMR was 
carried out over a wide range of pressure 0. The reso- 
nant frequency was found to be sharp in the low pressure 
phase but to broaden at higher pressure. Both the pres- 
ence of static magnetism and the broader distribution 
of frequencies (|M| distribution) is consistent with our 
proposal of BCCl for the high-pressure phase. 

Conclusions: A time reversal symmetry breaking he- 
lical spin crystal BCCl, disordered by weak impurities. 



is proposed for the high pressure 'partially ordered' state 
of MnSi. Both theoretical arguments from simple mod- 
els where this state naturally arises adjacent to the single 
spiral state, as well as favorable comparison with a vari- 
ety of experiments: e.g. neutron scattering and magnetic 
field studies lends support to this view. Future experi- 
ments that could probe the unusual properties of this 
phase e.g. magneto-conductivity and /j,SR are also pro- 
posed. Important puzzles that remain are the origin of 
NFL behavior and the lack of clear signatures for a finite 
temperature transition into the "partially ordered" state. 
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We have also checked the energies of some spin crystal 
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states outside of the six-mode model |l2||. The transition against all rivaling states tested so far. 

from the spin spiral to the BCCl state remains stable 



